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Longitudinal  Magnetoresistance  in  Polar  Semiconductors 
Robert  L.    Peterson 

The  magnetophonon  effect  due  to  resonant 
interactions  of  charge  carriers  with  optical 
phonons  in  nondegenerate  polar  semiconductors 
is  studied  when  several  elastic  scattering 
mechanisms  are  simultaneously  active  also. 
The  simple  model  of  the  band  structure,    and 
the  displaced  Maxwellian  distribution  function 
are  used.     Analytical  expressions  for  drift 
mobility  are  given  for  B  =   0  and  the  quantum 
limit,    and  are  compared  with  the  results  of 
other  techniques  where  possible.      The  numerical 
calculations  show  that  whenever  magnetophonon 
oscillations  are  seen  in  the  ohmic  regime,   maxima 
occur  near  the  resonance  magnetic  fields  for  all 
lattice  temperature  and  strengths  of  scattering 
mechanisms  which  compete  with  the  optical 
phonon  scattering.      This  is  in  contrast  with 
typical  experimental  findings,    in  which  minima 
usually  (but  not  always)  appear  near  the  resonance 
fields  for  ohmic,    longitudinal,    conditions.      The 
explanation  of  the  minima  must  therefore  be  sought 
in  the  use  of  a  distribution  function  in  which  carrier- 
carrier  scattering  is  less  important. 

Key  words:     Magnetophonon  effect;  magnetoresistance; 
semiconductors;  transport  theory. 

1.     Introduction 

This  Technical  Note  is  a  back-up  article  for  a  much  shorter 

paper  being  published  elsewhere.        We  here  give  a  rather  detailed 

exposition  of  the  calculations  only  indicated  in  that  paper,    and  provide 

some  further  discussion. 


2,  3 
In  earlier  articles,  we  developed  the  theory  and  showed 

the  results  of  numerical  computations  of  longitudinal  magnetoresistance 

for  the  simple  model  of  a  semiconductor,    using  the  displaced-Maxwell- 

ian  distribution  function  for  the  carriers,    and  scattering  of  carriers  by 

2 

optical  phonons  via  the  deformation  potential  ("nonpolar")     and  piezo- 

3 
electric  ("polar")     interactions.      Scattering  by  acoustic  phonons  via 

3 
the  deformation  potential  was  also  included  in  the  latter.        These 

computations   showed  oscillations  in  magnetoresistance,    electron  tem- 
perature,   and  acoustic  gain,    due  to  the  "magnetophonon  effect."     For 

the  longitudinal  ohmic  magnetoresistance,    such  oscillations  were  pre- 

4 
dieted  theoretically  in  1964,      and  have  since  been  seen  in  abundance 

5-16 
experimentally.  However,    the  experimental  observations  in  the 

longitudinal  ohmic  case  have  not  appeared  to  bear  out  the  theory  in 
more  than  a  cursory  way. 

The  experiments  typically  show  minima  in  the  ohmic  longi- 
tudinal magnetoresistance  near  the  resonance  fields,    i„  e.  ,    the  magnetic 
fields  for  which  the  longitudinal  optical  (LO)  phonon  energy  is  an  integral 
multiple  of  the  cyclotron  energy  of  the  carrier.     However,    in  certain 
materials  neither  the  maxima  nor  the  minima  lie  close  to  the 

resonance  points. 

4 
Gurevich  and  Firsov,      by  solving  the  Boltzmann  equation 

in  an  approximate  way  in  the  ohmic  region,    had  shown  that  for  pure 

polar  LO  scattering  the  extrema  of  the  longitudinal  magnetoresistance 

2 


oscillations  would  be  maxima  near  the  resonance  fields.      Our  com- 
putations based  on  the  displaced-Maxwellian  also  showed  this  for  both 

2  3  3 

the  nonpolar     and  polar     cases.      Gurevich  and  Firsov     further  deduced 

that  when  acoustic  phonon  scattering  competed  sufficiently  strongly 

with  the  LO  scattering,    these  maxima  would  convert  to  minima.      They 

evaluated  a  parameter     r    which  characterizes  the  relative  strength  of 

acoustic  to  optical  phonon  scattering.     When    T  is  greater  than  about 

4 
unity,    minima  should  be  observed.        Experiments  showing  minima 

have  been  reported  in  which    r  has  been  estimated  to  be  larger  than 

12 
unity,         supporting  the  Gurevich-Firsov  theory.      However,    in  other 

experiments  showing  minima,    it  appears  that     p    is  substantially  less 

7    8 
than  unity    '      (see  also  section  3). 

Further,    it  is  not  evident  just  how  the  maxima  for  pure 

optical  phonon  scattering  will  convert  into  minima  with  increasing 

strength  of  competing  scattering  mechanisms:     Do  the  magnetoresist- 

ance  extrema  shift  to  higher  or  lower  magnetic  field,    as  might  be 

10,  16 
inferred  from  the  more  polar  materials,  or  do  the  amplitudes 

of  the  extrema  simply  reverse  their  signs  with  relatively  little 
shifting  of  the  positions  of  the  extrema? 

These  unresolved  questions  have  prompted  us  to  begin  a 
more  detailed  theoretical  study  of  this  situation.     We  report  here  the 
results  of  an  investigation  in  which  the  displaced-Maxwellian  distribu- 
tion function  was  used  in  a  study  of  the  longitudinal  magnetore  si  stance 

3 


in  polar  semiconductors.      Four  scattering  mechanisms  were  considered 
simultaneously  in  varying  strengths:     the  polar  optical  phonon  scattering 
mechanism,    scattering  of  the  carriers  by  acoustic  phonons  via  both  the 
deformation  potential  and  piezoelectric  mechanisms,    and  also  scatter- 
ing due  to  screened  ionized  impurities.      Computer  studies  were  made 
at  several  lattice  temperatures.      We  have  used  the  simple  model  of  the 
band  structure,    that  is,    the  model  in  which  there  is  but  a  single  iso- 
tropic,   parabolic  conduction  band  centered  at  k=  0,    and  we  do  not  include 
Landau  level  broadening.      This  is  the  same  model  as  used  by  Gurevich 
and  Firsov. 

(If  we  had  included  nonparabolicity,    which  is  an  important 

requirement  for  the  much- studied  n-type  InSb,    the  displaced- Maxwellian 

17 
could  not  have  been  used.      The  same  techniques  used  in  deriving       the 

displaced-Maxwellian  for  a  parabolic  band  result  in  a  very  much  more 
complicated  distribution  function  in  the  nonparabolic  case.  ) 

Our  reason  for  choosing  this  particular  distribution  function 
is  that  it  is  the  only  one  which  at  present  seems  to  be  amenable  to 
extensive  numerical  computation  for  polar  semiconductors.      Use  of 
the  displaced-Maxwellian  of  course  implies  a  good  deal  of  carrier- 
carrier  scattering  which  establishes  the  distribution.      Conditions  for 

1  8 

the  validity  of  this  distribution  are  discussed  by  Conwell,         who 

references  earlier  work.      Reference  2  contains  further  references  to 


work  based  on  the  displaced-Maxwellian.      If  it  is  felt  that  carrier- 
carrier  scattering  is  not  predominant  in  a  given  material,    then  of 
course  the  displaced-Maxwellian  results  should  not  be  taken  too 
literally. 

In  the  next  section,    we  establish  the  basic  formulas 
corresponding  to  the  model  chosen.      We  derive  formulas  for  ohmic 
drift  mobilities  for  the  various   scattering  mechanisms  for  both  B  =   0 
and  the  quantum  limit,    that  is,    the  limit  in  which  only  the  lowest 
Landau  level  is  considered  to  be  occupied.      These  formulas  are 
compared  with  ones  obtained  by  solution  of  the  Boltzmann  equation, 
where  possible.     In  section  III  we  show  several  graphs  of  ohmic 
magnetoresistance  at  various  scattering  strengths  and  temperatures. 
We  also  exhibit  some  hot- electron  magnetoresistance  curves.      Dis- 
cussions of  the  results  are  given  and  some  conclusions  are  drawn. 

2.       Theoretical  Basis 

The  displaced-Maxwellian  distribution  function  contains  two 

parameters,    which  are  the  carrier  drift  velocity    v       and  the  carrier 

temperature     T     =   l/kg    .      These  are  determined  self- consistently 
e  e 

through  two  coupled  equations  which  describe  the  constancy  of  the 

2 

carrier  momentum  and  energy  in  the  steady  state,      or  in  other  words, 

which  balance  the  drift  and  scattering  terms  of  the  first  and  second 
moments  of  the  Boltzmann  equation. 


The  scattering  interactions  for  the  three  carrier-phonon 
mechanisms  considered  here  may  all  be  written  in  the  form 

Eic±'L 
"q     b    e  +  Hermitian  conjugate,  (1) 

where    b       is  the  phonon  destruction  operator  for  mode    q,     and   r     is 
q  «•  «• 

the  carrier  position.      The  index    a   indicates  a  specific  scattering 
mechanism,    and  polarization  indices  have  been  suppressed.      The 
summation  is  over  all  phonon  wave  vectors    q  in  one  or  more  polariza- 
tion  branches.      The  coupling  coefficient    q      is  generally  anisotropic, 
and  it  is  customary  to  substitute  an  averaged  quantity  for  it. 

The  scattering  interaction  between  a  carrier  and   N 
randomly  distributed  screened  ionized  impurities  may  be  written 


I 

19 


iq  •  (r-R) 
p    ^-5 5—     '  W 

lmP  2-2        q2+r    ~2 


where,    in  SI  (or  MKS)  units, 


b=  e2/xe„      .  (3) 


o 


Here    K     is  the  dielectric  constant  and     e        is  the  permittivity  constant. 

o 

The  summation  in  eq  (2)  is  over  all  impurity  positions    R   and  Fourier 

components    q.      The  screening  length  is    r   .     An  obvious  extension  of 

19 
the  work  of  Argyres  and  Adams        to  the  hot- electron  regime  yields 

the  expression 


2  -2 

q        =    r  =  be  n     ,  (4) 

s  s  e    c 


where    n      is  the  carrier  concentration. 
c 

Subsequent  expressions  are  written  for  electrons,    of  charge 
-e    with    e   >  0.      Drift  mobility  is  defined  by 


v   (E,  B)=   -n  (E,  B)  E.  (5) 

d 


The  magnetoresistance  is 

n   (E,  0)u(E,  0) 


A  p      _     P(E,  B)    _j  _ 


PQ  P(E,  0)  nc(E,  B)  U(E,  B) 


(6) 


In  this  paper  we  shall  ignore  any   E    or    B    dependence  of  the  carrier 

concentration  n   ,    so  that  under  ohmic  conditions,    the  magnetoresist- 
c 

ance  is  calculated  from  the  ratio  of  drift  mobilities.      For  nonohmic 
constant- voltage  conditions,    the  mobility  ratio  in  eq  (6)  is  calculated 
from  the  ratio  of  drift  velocities,    while  for  nonohmic  constant- current 
(or  constant  drift  velocity)  the  mobility  ratio  is  given  by  E(B)/E(0), 
where  E(B)    is  the  B-dependent  electric  field  across  the  sample. 

We  consider  the  B  =   0  and  B  /  0  cases  separately,    treating 
the  zero  magnetic  field  case  first. 

2..  1     B  =   0 

The  complete  Hamiltonian,    with   E    in  the  z- direction,    is 

H  =  p   /2m'"  +  HT    +   13- V    +  eEz  ,  (7) 

L.         a    a 


where   p    is  the  electron  momentum,    m     is  its  effective  mass,    H      is  the 
***  L 

lattice  Hamiltonian,  and  the  interactions  V     have  been  defined  above.   The 

a 


displaced- Maxwellian,    normalized  to  unity  is 

3/2 
- 1  2        *  2       —  2  * 

p(k)  =  Q     (2nBh/m    )       exp  [  -  8  h  (k-k)  /2m    ]    ,  (8) 


where    Q   is  the  volume,    and    hk=m   v   .      The  steady- state  equations 

«*  "*d 

2 
for  electron  momentum  and  energy  are 


momentum:     eE  =  4L1  S  ,  (9) 

a    m,  a 

energy:  eEv,  =    Zj    S  .  (10) 

d        a      e»  a. 

The  left  sides  of  these  equations  come  from  the  drift  term  of  the 
Boltzmann  equation,    and  the  right  sides  are  the  scattering  terms.      The 
summations  are  over  all  scattering  mechanisms,    and  we  shall  not 
consider  possible  correlations  between  the  various  scattering  processes. 
Also,    we  take  the  collisions  with  impurities  to  be  perfectly  elastic,    so 

that  the  impurity  scattering  makes  no  contribution  in  eq  (10).      The 

2 
scattering  terms  may  be  written 


Se.a=(h2/2m*»      kEk'<k'2-k2»Wk'akP<ii'-  (12) 


The  transition  probability  per  unit  time  for  electron- phonon 
scattering  is 

a 

.  =   2^     P  (\  r_ 

q  a  q      <*,  q 


Wk'k=   ^     PHnJ)wa(|n^}    k',     |njk)      .  (13) 


The  summation  here  is  over  all  phonon  numbers    n       and  n  '     in  the 

q  q 

sets     fn    1   and     f  n'    }.      The  lattice  distribution  function     n({n    })  yields 

•  q  q "  q 

the  Planck  factor  for  the  number   n  of  thermal  phonons  in  mode     q 

q»  a  «. 

and  branch    n : 


n         =   [exp(B   h«)       )-l]  ,  (14) 

q,  a  o      q,  a 

where     uu  is  the  frequency  of  phonons  of  type    &,    and  T    =  l/kR     is 

q»  a  L  o 

the  lattiqe  temperature.     Also, 


w    -(2TT/h)|  <{n'  }  k'l  V     |{n}k>|  X 

q  -       a.       q  *- 


6   [£(ln'  jk')-M{n    1  k)  ]    .  (15) 

q    —  q   «rf 


The  energy    g,   is  given  by 


2    2 

£({n    }k)=  ^4r    +    E  h^        (n     +    I).  (16) 

q    ~  2m-  q  q,a      q         2' 


For  ionized  impurity  scattering,    the  transition  probability  is 
simply 


W1   /limp=(2TT/^)|<k/|V.         I  k>   |  2   6  [h2(k2-k'2)/2m*   ]  .         (17) 

iA.       fa. 

The  momentum  and  energy  scattering  terms  for  carrier- 
phonon  processes  reduce  to 

1 


I2 

m,  a  *'0  a  '      "»j  a 


S  =  A 


j      dq  q      I  q      I        J     dx  xF    (q,x),  (18) 


*°  2 

S         =   A     f    dq  q     uu  |   q       |  f    dx  F    (q,  x),  (19) 

e,  a  «L  q»  a       a  J  a 


after  integrating  over  carrier  coordinates.      Here 


A=  (2n3  m')2  Q/4TT     h   ,  (20) 


F      (q,x)=(n         +  1)  e  _Tla+  -Ti      ^e^^'     ,  (21) 

a  q»  a  q, a 


*2 

n         r  m    (jo 

=   -hr-     5-  ±  :   4'~    Tkx|      .        (22) 

2m' 


In  these  equations,     x  =   cos    9    where     9  is  the  angle  between   q  and  E. 
If  elasticity  is  assumed  for  the  acoustic  phonon  processes,    the  term 

m    iw         /hq  of  eq  (22)  does  not  appear.      However,    in  the  numerical 

q,  a 

work  the  term  causes  no  real  complication,    and  so  we  have  not 
assumed  elasticity  there.      Neither  have  we  assumed  equipartition, 
but  have  used  the  full  Planck  factor,    eq  (19),    in  the  numerical  work. 
(Equipartition  becomes  increasingly  inaccurate  with  increasing    E.  ) 


10 


We  now  consider  the  four  mechanisms  separately  in  more  detail. 

2..  1.  1   Deformation  potential,    acoustic  phonons. 

20 
The  coupling  coefficient  can  be  written 

^d,acl2=Elh<l2/2Pn%'  (23) 

where    E      is  the  deformation  potential  energy  and    p  is  the  mass  density. 
Only  longitudinally  polarized  acoustic  phonons  are  involved,    and  we 
use 

0Jq=U^q,  (24) 

21 

where    u       is  the  longitudinal  speed  of  sound  averaged  over  directions. 

The  E  =   B  =   0   drift  mobility  expression  is  readily  found  to  be 


1 
2    4  2 

3eputh       /  Zt?      \ 


16  E  \m      k    T         / 

1  J_i 

This  is  obtained  from  the  ohmic  limit  of  eqs  (9)  and  (18),    using  elasti- 

22 
city  and  equipartition,    and  has  been  derived  previously  by  Sladek 

23 

and  Paranjape  and  deAlba.  The  mobility  as  given  by  eq  (25)  is 

smaller  than  that  obtained  by  solving  the  Boltzmann  equation  for  de- 

24 
formation  potential  acoustic  phonon  scattering  alone,        by  the  factor 

9tt/32.      The  reason  is  that  the  displace d-Maxwellian  implies  carrier- 
carrier  scattering,    which  further  limits  the  mobility. 


11 


2.1.2    Piezoelectric,    acoustic  phonons. 


For  a  given  direction  of  propagation    q,    the  contribution  of  a 


given  polarization  branch  to    |q      _^|      may  be  written 


2  .  ...       25 

ac1 


K2bhu)    /2nq2,  (26) 

q 

where   K    is  the  dimensionless  electromechanical  coupling,    and  b    is 
given  by  eq  (3).      Both  quantities  are  directional  dependent.      We  use 
an  average  value  for  the  dielectric  constant  appearing  in  b.      The  ex- 
pression (26)  is  summed  over  longitudinal  and  transverse  modes,    and 

21  26 

averaged  over  directions        in  the  manner  of  Meijer  and  Polder.  We 

indicate  the  result  symbolically  here  by 


S,ac|Z=Bp/Q"  <27' 


and  give  additional  details  regarding  computation  in  section  III. 
The  E  =   B  =   0  drift  mobility  is  found  to  be 

3eu      h  /  3       \ 

av  /        2tt  \ 

"o  =  ^i [—3 )  (28) 

p  \m      kTL / 

22 

in  the  elastic  equipartition  limit,    a  result  given  ea.rlier  by  Sladek. 

As  in  the  preceding  case,    the  result  (28)  differs  from  the  Boltzmann 
equation  result  only  in  being  smaller  by  the  factor  9tt/32. 
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Here  we  use 


2.1.3    Polar,    optical  phonons. 
27 


1 2  2  *    2 

q  =   2Tfh    eE   /Qm     q     .  (29) 

^op  o  n  v 


E     has  units  of  electric  field  and  is  given  by 
o 


m    eW 


euJ  / 

s  h      W         x     y  ' 

O  v      ^°  O  / 


E     =^ —       —    -     —       »  (30) 

o      4rre 


where     K      and    k      are,    respectively,    the  (directionally  averaged)  di- 

o  co 

electric  constants  for  zero  and  infinite  frequencies,    and    ^       is  the 
longitudinal  optical  mode  (LO)  frequency  at    k  =    0.      We  make  the  usual 
assumption  that  the  LO  frequencies  at  all    k   are  equal  to    uj    . 

The  phonon  scattering  expressions  (18)  and  (19)  are  not  the 
simplest  to  use  for  optical  phonon  scattering.      It  is  easier  to  integrate 
first  on  the  phonon  modes    q,     then  over  the  angular  parts  of  k,    re- 
suiting  in 

2     .oc  2 


S 
m 


2eE    n/\2_Yv       C  _yx      I       v        Y 

o    o/y\  1  J      Yo  -  ' 

,op= v~  S^J      e  JodXe  e  ^     [xs+tn(x+s)] 

(31) 
X    C(2yvs)    +    —     [xs  -   ta(x+s)]C(2Yvx)k 
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-  Y  x 


/2h^oy    \  -W  « 

'     I  — ^ J  n   e  I     ( 

o\      *  /         o  J 

^m     tt       /  o 


S  =  4eE     I  ; )      n   e  I     dx  e  x£n(x+s)  X 

e ,  op  o  \       ->: 

^m     tt 


rYo~Y  i 

|e  sinh(2Yvs)  sinh(2Yv  x)  J 

••  2YvS  "         2vvx         J' 


2        \ 
where    s  =  (x  +1)    ,    and 


The  E=  B=  0  drift  mobility  is  found  to  be 


,  2  v  / 

3      /         ^o      \  sinh  (V  2) 


(32) 


C(y )  =   cosh  y  -  y      sinh    y  ,  (33) 


*    2  2 

v    =   [m  vd/2ha)Q]      ,  (34) 


Y 
\    =   fe    ""-I]"1   ,  (35) 


Y      =   p  hi«  ,      Y=   3  huo  (36) 

o  o       o  e       o 


/  TThuuo 

^O-Wi— ~ j)         kTVz)^    '  (37) 

o   \  2m    Y      /  1    Y  o7 

o 

28 
first  derived  by  Stratton.  Properties  of  the  modified  Bessel  function 

K     may  be  found  in  Ref.    29.      The  high  temperature  form  of  eq  (37)  is 

9TT/32    times   smaller  than  that  deduced  from  the  Boltzmann  equation, 

as  in  the  acoustic  cases  (equipartition  assumed  there  is  a  "high" 
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temperature  approximation).     At  lower  temperatures  (Y    >  0.  36)    it  is 
well  known  that  the  Boltzmann  equation  cannot  be  solved  directly  for 

the  polar  optical  interaction  because  of  the  breakdown  of  the  relaxation 

30 
time  approach.     Howarth  and  Sondheimer       used  a  variational  technique, 

which  when  corrected  by  the  effective  ionic  charge  in  the  manner  of 

31 
Ehrenreich,         gives  precisely  eq  (37)  in  lowest  order  approximation. 

2.1.4    Screened  ionized  impurities. 

From  eqs  (2)  -  (4),    (9),    and  (17),    one  can  write 


S 

m,  imp 


2       ,„.       .1 

L2  2 


-n    b  /  2B      \ 

— —  I — 2.  I 

16  v,  13     */ 

d  ^rr  m  / 


(38) 


1 


1 
dx    x""1   [£n(l+x)-x/(l+x)]e  ~5"C[2r(gx)2]  , 


where    n    =  N  /q     is  the  impurity  concentration, 

g    =   3    e    /4  ,  (39) 


e     s 


e       =    h   q  2/2rrT=h  8  n   b/2m"~,  (40) 


e    c 


r2  =     8    m*v   2/2    ,  (41) 

e  d 

and  C(y)  is  given  by  eq  (33).     In  the  ohmic  region,    the  integral  in  eq  (38) 

32 

can  be  expressed  in  terms  of  the  exponential  integral         E ..(g): 
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Vs)  -  J" 


•  gx        -  "• 


E1(g)=     i!      dx  e  x  .  (42) 

The  E=B=0  drift  mobility  is 


3.  3_     3 


U 


/2tt  k   T        \        r 

o-^i—i^-)    K> 

nb         x         m  /        u 


e       °    E^g   )-I 


,(43) 


where    g       is  the  E   -»   0   form  of   g.     Equation  (43)  is  apparently  new. 


o 


19 

Solution  of  the  Boltzmann  equation       in  the  ohmic  limit  yields  a  some- 
what different  integral  than  that  to  which  eq  (38)  reduces,    which  cannot 
be  expressed  in  terms  of  exponential  integrals.     However,    for   g     <<:1, 
which  is  typical  in  nondegenerate  semiconductors,    eq  (43)  may  be  sim- 
plified by  using  E    (g    )  «    In  (l/g    ).      Using   g     <<    1    in  the  Boltzmann 

19 
equation  result,         one  obtains  the  same  expression  as  that  to  which  eq 

(43)  reduces,    except  that  the  drift  mobility  in  the  displaced  Maxwellian 

formulation  is  3rr/32  times  smaller. 

2.2    Finite  Magnetic  Field 
In  the  familiar  "Landau  gauge,  "  the  Hamiltonian  for  E    and    B 
both  in  the  z- direction  may  be  written 

2  2  *  2  *  V* 

H  =    [p        +  P      +(p   +m    uj     x)    ]/2m    +  H  _  +  2-jV   +eEz    ,  (44) 

X  Z  y  C  J_j  rv        Ct 

where 

0)     =   eB/m  (45) 

c 


lb 


is  the  cyclotron  frequency.      The  carrier  energy  in  the  plane  wave 

states    k   ,    k   ,     and  Landau  level   n  =   0,    1,    2,    ...    is  thus 
y      z 

£     (k   )  =  h2k   2/2m'":  +  huu    (n+  i)  .  (46) 

n     z  z  c  2 


We  have  associated  the  energy  degeneracy  with  the  y- direction,    with 
N     degrees  of  freedom.      The  displaced-Maxwellian  may  be  written  as 

y 

the  product    p       n  (k   ,  k    ),      where 
n  y      z 

1 

nf*      V    1         2^       /BehY         ^'e^1^^)2]    ,(47) 
y    z  \2rtm    /  2m 


=  (l-exp(-B    h^   ))     exp(-nghu)    )    .  (48) 

e      c  e        c 


The  distribution  function  is  normalized  such  that 


CO 

£         E      (       dk      p     p(k,k)    =    1, 

1=n      k     J-„      z    n       y     z 


(49) 


n=0       k 

y 


and  L      is  the  length  of  the  sample  in  the  z- direction. 

in  •    i» 
The  matrix  elements  of  e    A     «•%    needed  for  both  phonon  and 

impurity  scattering,    are 


iq«  r 

'tiK  k  >=  o, 
v    z  I 

y         y     y-        z         z       z 


<n'k    'k    'le  Ink  k  >=.fi.      ,   ,  6  .      /  .  J    ,    (q  ,k  +q  ,k   ),    (50) 

y      z    I  '      y    z  k^'-k^q        k^  -k^,  q^    n   n  nx      y      y      y 
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where 


lq       x 
x 


J     /„(q   ,  k   +q   ,  k    )    =   f      dx  e  u       /  x+  k  +q     \u     (x  +  k  ),      (51) 

nnxyyy         «'  n. ' !  y      y  J    n  \  y         / 


m    uu   /h  m  'uo  /h 

c  c 


and  u   (z)  is  the  harmonic  oscillator  wave  function 
n 

_i 
n  *  4  2 

i   H   (z)        /m    (i)    \  -z  /2 


i  H  (z)         /m    uo    \  -z  /2 


(2     n!) 


and  H   (z)  is  the  Hermite  polynominal  of  order   n.     Now   J     ,     is  a 
n  n  'n 


function  only  of   \  ,      where 


2  * 


2  2 

and   q         =   q     +q_     „      One  finds 


\ 


I  J     ,     ]2=   e"Xf    ,(X)     ,  (54) 

n    n  n    n 


where 

r  n  s       "1 2 

f    ,    (X)=  n!nMXn'-n    I    E  ,,["V   -    .,  .(nfen).      (55) 

n'n  L  s=o        s !  (n- s)!  (n   -n+s)!    J 

If  equipartition  and  elastic  collisions  are  assumed  for  the  de- 
formation potential  acoustic  phonon  scattering,    the  greatly  simplifying 
result 
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CO  00  CO 

h(2flm*tti  )_1     J        dq         J         dq    |J        |     =    f   dXe    \        (X)  =1  (56) 
''-od        x       J-o=  y       n'n1  J  0  n'n 

2 

can  be  used,    as  it  can  be  with  the  nonpolar  optical  phonon  scattering. 

However,    for  the  other  scattering  mechanisms  considered  here,    this 
result  cannot  be  used,    and  so,    as  with  B  =   0,    we  have  not  assumed 
equipartition  or  elastic  collisions  for  the  acoustic  phonons,    as  no 
significant  computational  time-saving  is  thereby  achieved. 

The  steady- state  equations  for  carrier  momentum  and  energy 
can  be  written  in  the  forms  eqs  (9)  and  (10),    respectively.     For 
carrier-phonon  interactions,    the  scattering  terms  are 


1  2 


Sm,a  <B>  =   "A   J  dVXl         J    dq|>qjLl^a  |      \<V  q,). 


(57) 


CO  -  1  oo 

Seja(B),-A      j       dqjqj  /     dqi(«U|-iIo|2%aRa(qiifqB).    (58) 


where     A    is  given  by  eq  (20),     and 


R    (q    ,  q   )=         ^n  p      e       f    ,   (X) 

a    U    nz      n,  n  '=  0    n  '  n'n 


i  i 

(n  +1)  e      a  +  -n  e  j  ,     (59) 

q,  a  q,  a 


.  R    h        I    q  m'uJ(n-n')        m   (l)  J     2 

=  l£_l^+ £ ±      ,  q'a    T  *j        . 

a  i      o    *         2  hq„  i^q. 


2m  hqz 


(60) 
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For  screened  ionized  impurities,    the  momentum  scattering 

term  may  be  written 

,        2 

2   -  r         ..  °°  °° 

-n  bhui  Be  f  2       f  _ 

|B>'      ~   3/26 J0dye  J„dXe       [VU'yl 


s 

m,  unp 

4tt 


(61) 


+  Rimp(^-y)J  ' 


where 


co  n 

R. (\,y)=      2-  p     |      E     f  ,(l)c      (X,  y)  sinh  (2ry) 

-1=0 

(62) 


imp  n=  0    n  L    „_  n      n,  n-£ 


+      Sfn)n^(X)exp(-,,Reh-c)c^(X,y)sinh(2ry^J   , 


y      =        (y    Uphffl    )2,  (63) 

■C  e         c 


c 
I 


..      -  2 

(\,y)  =  y^'V^  y)[DU)  +  (y^+y)2]       .     (64) 


D    (  X)  =   8    h«J  X  +  4g     ,  (65) 

e      c 


and   g   and   r    are  given  in  eqs  (39)  and  (41),    respectively. 

We  now  look  at  some  results  in  the  quantum  limit,    which  we 
define  here  to  mean  that  the  term     e    e         c    is  negligible  compared  to 
unity,    so  that  all  carriers  are  in  the  lowest  (n=  0)  Landau  level. 
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Z  .  2.  1     Deformation  potential,    acoustic  phonons. 

The  two  acoustic  phonon  mechanisms  can  be  treated  together 

i—    1 2  Pa     , 

up  to  a  point,    by  writing     |  q  =   B    q         /Q,    where    p     =  +1     for  the 

a  a  a 

deformation  potential  interaction,    and    p    =  -1    for  the  piezoelectric 

a 

interaction.      The  momentum  scattering  term  (57)  may  then  be  written, 
in  the  equipartition,    elastic,    ohmic,    and  quantum  limits. 


:<  co  p  -i 

-8AB    m      v    ou      f       ,       -y  f  -X     a 


Smj    (QL)  =  a d_c  dy  e        V       dXe 

J  _.  ^  o  •»  o 


(66) 


1  e0ua  c 


and 


=  [smV  P  Qh  2  +  2m^c X/h  J  2     .  ( 67 ) 


The  ohmic  drift  mobility  for  deformation  potential  scattering  is  then 

3,3  2 

rr    eh     pu. 

^QL=  2        ~      *5  i  <68) 

Ej    UJc(2m     kTL) 


Comparing  with  the  B  =   0  case,    eq     (25),    we  see  that 

U>QL=3hu)c/8kTL    •  (69) 

Apart  from  a  possible  variation  of  carrier  concentration  with  magnetic 
field,   the  longitudinal  magnetoresistance,    eq     (6)     thus  increases 
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linearly  with  field  in  the  ohmic,    quantum  limit.      Equation  (69)  agrees 

19 
with  the  Boltzmann  equation  result  without  carrier- carrier  scattering, 

except  that  a   9    replaces  the    8    in  eq  (69). 

2..  2.  2    Piezoelectric,    acoustic  phonons. 

32 
Using   p      =   -1,    and 

a 


i 


dx  e       X  E   (x)  =  a  "      £n(l  *   a)  (70) 


in  eqs  (66)  and  (67),    one  finds 


1 

("3tM'c-i 


4h    eu 

av      I  Li        c-1 


where 


UQL  B    uu  \_     *3      J        In  c  (71) 

p    c  Zm 


c  =    3    hiu  /4  .  (72) 


Comparing  with  eq  (28),    we  have 


^o  3 


QL 


c  In    c 
2        c-l 


(73) 


Thus  the  ohmic  longitudinal  magnetoresistance  increases  logarithmically 
with   B    at  sufficiently  high  field,    apart  from  any  changes  in  carrier 
concentration.     Except  for  the  numerical  factor,    eq  (73)  agrees  with 
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that  obtained  from  the  Boltzmann  equation  without  carrier- carrier 

33 
scattering,  for    c  >>    1. 


2  .  2.  3    Polar,    optical  phonons. 

In  the  ohmic  quantum  limit,    the  momentum  scattering  term  (57) 

becomes 

l 


Li 


where 


I(TL,  B)  =  J  Q  dx  E^x)  exp  [-(c-1)  x- yq2/(16  ex)]  ,         (75) 


and   c   is  given  in  eq  (72).      The  ohmic  drift  mobility  is  thus 

4  sinh  (Yo/2) 


l, 
3„     3  .      2 


WQL"V    ,      *      /  Et„     K'i'.B)        •  (?6) 

2m        /  o       c  Li 


For   c    sufficiently  greater  than  unity,    it  is  evident  that  the  term 

? 
exp  ["-   V      /(16cx)l    in  I(T    ,  B)  makes  no  essential  contribution  to 

L       '  o  L 

I(T    ,  B)  apart  from  a  temperature- dependent  multiplicative  factor. 


We  thus  write,    using  eq  (7  0), 


I(T_,B)«  r(YJ  (c-1)"1    in    c    .  (77) 
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That  eq  (77)  is  a  good  approximation  can  be  seen  from  figure  1,    where 
the  exact  integral,    evaluated  numerically,    and  eq  (77)  are  plotted 
against   c    for  three  temperatures.      The  factor  r(y    )    is  there  taken 
to  be  exp  (-0.44  y    ),    chosen  strictly  empirically.      Thus,    to  a  good 
approximation,    we  can  write  the  drift  mobility  as 

L,   \      c-1  sinh  (y    /2) 

—}    —, —  o  76  a) 

-    /    c  -In  c      — : r—      »  v  ' 

2m      v  E     r  (y   ) 

o  o 

and  the  ratio  of  mobilities  as,    with,  eq  (37), 


o  3r(Y   )  cine  (78) 


o 


UQL  y^^yJZ)       c-1 


The  ohmic  longitudinal  magnetoresistance  thus  varies  approximately 
logarithmically  at  sufficiently  large  B,    apart  from  any  carrier  con- 
centration variation. 

2  .  2.  4    Screened  ionized  impurities. 
The  momentum  scattering  term  (61)  reduces  to 


2 

2    -  r  f  2 

S  (QL)="nibe  |        dy  e"Y    sinh(2ry)z"1eZE ,(z)  ,      (?9) 

m,  imp  -r- J  2 

0     5/£  o 


in  the  quantum  limit,    where 
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c  =  Xo^c  /4aJo 


Figure  1.      Comparison  of  the  integral  I(  T     ,  B)  defined  in  eq  (75' 
(solid  line)  to  the  approximate  form  given  in  eq  (77) 
(dashed  line). 


Z5 


z  =  4(g+y    )/b   h<«      , 


(80) 


r   is  given  by  eq  (41),    and  E    (z)  is  the  second-order  exponential 
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integral.  In  the  ohmic  limit,    the  drift  mobility  is  found  to  be 

3 


QL 


J^~e 


8\2    eg   (TTkT     ) 

O  J_i 

nIb2m^[F(go)-F(go/c)] 


c-1 


(81) 


where 


F(x)  =     x   e      E^x)    . 


(82) 


With  eq  (43),    the  ratio  of  B  =   0  to  quantum  limit  mobilities  is 


o 


3c  F(gQ)  -  F(gQ/c) 


UQL  2(c-l)      (l  +  go)  F(go)-go 


(83) 


The  longitudinal  magnetoresistance  due  to  ionized  impurities  thus 

saturates  at  sufficiently  high  magnetic  field.      The  Boltzmann  equation 

19 

result  of  Argyres  and  Adams,         which  omits  carrier-carrier  inter- 
action,   also  shows  this,    but  with  a  different  saturation  value.      For 

example,    if    g      <<  1,      and    c     >  1,    then  using  E  (x)    «   -y     -£n  x 
o  ~  IE 

where      y      =   0.5772  ...    is  Euler's  constant,    one  obtains  from 
eq(83) 
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M  0        3        y  E  +  to  g0  +  (c-lf    in    c 


1Q 

to  a  good  approximation.      The  result  of  Argyres  and  Adams        reduces 


to  2[tn(3/go)-l]_1(l  +  3/c)"1, 


3.       Numerical  Results 


In  this  section  we  exhibit  some  results  of  our  numerical  v-ork. 

For  this,    it  is  necessary  of  course  to  choose  values  for  the  parameters 

such  as  m    ,    oj    ,    E,,    etc.      We  have  elected  to  use  those  which  seem  to 
o        1 

be  representative  of  n-InSb,    which  is  probably  the  most  extensively 
studied  material  to  date  regarding  magnetophonon  effects.      However, 
it  should  be  kept  in  mind  that  InSb  has  a  very  nonparabolic  central 
conduction  band,    whereas  our  model  assumes  a  parabolic  conduction 
band.      Further,    the  displaced-Maxwellian  distribution  function  implies 
a  degree  of  carrier- carrier  scattering  which  probably  exceeds  that 
in  most  samples  studied  so  far.      Thus,    the  numerical  work  presented 
here  should  not  be  understood  as  an  attempt  to  describe  n-InSb.     In- 
stead,   this  study  was  undertaken  to  determine  the  extent  to  which 
temperature  and  elastic  scattering  could  shift  the  magnetophonon 
structure,    within  the  displaced-Maxwellian  formulation. 
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Figures  2,    3,    and  4  show  curves  of  ohmic  magnetoresistan.ee 
versus  magnetic  field  (expressed  as  oj    /iw   )  for  three  temperatures 
corresponding  to    y      =   Q/T      =  h  uj   /kT      =   14.0,    3.65,    and  2.  60, 

O  J_j  O  L-i 

respectively.      With    9  =   281  K,         these  values  correspond  to  lattice 
temperatures  of  20  K,    77  K,    and  108  K.      We  have  also  used  E     =  470 

V/cm,    and  m     =   0.  0159  m    ,    which  is  the  value  of  the  so-called 

'  o 

11 

magnetophonon  mass  in  n-InSb.  The  value  of  the  deformation 

34 
potential  E    has  been  uncertain,    although  Tsidilkovskii  and  Demchuk 

have  recently  made  a  strong  argument  for  a  value  of  about  3  0  eV. 

They  reference  much  previous  work,    but  omit  the  work  of  Nill  and 

McWhorter        (4.  5  ev),    Whalen  and  Westgate        (16.  2  eV),    and  Martin 

37 
and  Mead        (16  eV  or  higher).      In  the  figures  we  have  used  several 

values  for  E    .      Also,    we  use  the  Nill  and  McWhorter  value  of  0.  06 

coul/m     for  the  piezoelectric  constant  e      ,    which  is  close  to  that  of 

Whalen  and  Westgate.      The  piezoelectric  coupling  coefficient  for 

acoustic  phonon  scattering  is  obtained  in  the  manner  of  Meijer  and 

2b 
Polder       by  considering  all  polarization  modes  and  averaging  over 

directions.      The  result  is 


B     =   2b2e142hQ/pe2    ,  (84) 


where  B     appears  in  eq  (27).      The  averaged  quantity   Q   is 
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Figure  2.       Ohmic  magnetoresistance  at   y      =   14.  0  (20  K  for  InSb). 

Values  of  E     (eV),    e        (coul/m    ),    impurity  concentration 

-3 
(cm      ),    and   T    are,    respectively:   (a)  0,    0,    0,    0;  (b)  4.5, 

0.  0b,    4x10    3,    0.  27;  (c)  lb.  2,    0,    0,    0.  01;  (d)  lb.  2,    0, 

15 
10      ,    1.5b. 
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Figure  3.       Ohmic  magnetoresistance  at   y      =   3.  65  (77  K  for  InSb). 

Values  of  E     (eV),    e        (coul/m     ),    impurity  concentration 

-3 
(cm      ),    and    T    are,    respectively:  (a)  0,    0,    0,    0,  ;  (b)  4.  5, 

0.  06,    4   x  1013,    0.  04;  (c)  16.2,    0,    1015,    0.43;  (d)  30,    0.  0b, 

1015,    0.  58;  (e)  0.  58;  (e)  30,    0.  06,    4   x  1Q13,    0.  23. 
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Figure  4.       Ohmic  magnetoresistance  at   y     =   2.  bO  (108  K  for  InSb). 

Values  of  E     (eV),    e        (coul/m     ),    impurity  concentration 

_  3 
(cm      ),    and    T  are,    respectively:   (a)  0,    0,    0,    0;  (b)  4.  5, 

13  15 

0.0b,    4x10      ,    0.  03;  (c)  lb.  2,    0,    10      ,    0.  31;(d)30,    0.0b, 

10      ,    0.51;  (e)  30,    0.0b,    4x10      ,    0.30. 
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(3ul 


4u~ I   /13     '  (85) 


where    u    and  u       are  longitudinal  and  transverse  sound  velocities  for 
propagation  in  the  (1,    1,    1)  and  (1,    1,     0)  directions,    respectively, 
given  by 


pUl2  =   I   (C11  +  2C12  +  4C44  +  4ei42^£o)  '  (86) 


PU22   =     C44+ei42/Keo    •  (8?) 


The  elastic  constants  are  given  the  values     c,  ,   =   6.  66,    c        =   3.  35, 

B  11  12 

10        ,2 
and   c=   3. 14x10       nt/m    .      We  also  chose     K  =   17.8  and    D=  5.  78 

44  '  M 

3 
gm/cm    .      The  calculated  values  of  u     and   u      are  3900  and  2300  m/sec, 

respectively,    and  the  grand  average   u        over  directions  and  polariza- 
tions is  2800  m/sec,    which  is  used  in   oj     =  u      q.      For  the  deformation 

q         av 

potential  scattering  we  used   u    =   3700  m/sec,    which  is  the  average 
longitudinal  sound  velocity. 

In  evaluating  the  screening  terms    g   and   g    ,    we  have  set  the 
carrier  concentration  equal  to  the  impurity  concentration.      Finally,    we 
have  summed  over  as  many  Landau  levels  as  necessary  to  achieve 
about  0.  1%  accuracy.      In  some  cases,    as  many  as  14  levels  were 
included. 
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4 
We  have  calculated  the  Gurevich-Firsov       T  (see  the  discussion 

in  the  Introduction)  for  each  of  the  curves  in  figures  2-4.      Extending 

their  result  to  include  all  scattering  mechanisms  which  compete  with 

the  carrier- optical  phonon  scattering,    and  assuming  that  the  respective 

resistivities  are  simply  additive,    we  obtain 


(      h      )  *E 

I  IT  uu    kTT        /  a 

\         o      L     / 


r=     T    ^c      n  UJ    kT         /         ~    Pa 


The  three  dimensionless  resistivities     p       in  the  present  case  are 

a  c 


?  L 

2      *  /  3         3     w   2 

16E     m  /  ha;    k    T  \ 

(def.pot.,ac.)=    2 5~  V °   5  /      >  (89> 

3eE   tl     pu  V  TT 

o       M    £ 


a      **  2         2- 

8m 


'b    e         Q            /h«J    kT      \  2 
p(piezo.  ,ac.  )=   — I  V/       '         (90) 
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K2  /         ,  * 

n  _b  /       h  u) 


n  b  /       hu)         \  &o 
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F- values  were  calculated  at  the  first  resonance  point,    that  is  at  w  =w 

c       o 

As  is   seen  from  the  figure  captions,      F    is  typically  below  unity  even 
in  the  most  favorable  cases.      Thus  according  to  Gurevich-Firsov 
theory,    resonance  point  maxima,    or  perhaps   simply  very  small 
amplitudes  of  oscillation,    are  to  be  expected  in  n-InSb.     However, 
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as  pointed  out  in  the  introduction,    magnetophonon  oscillations  in  n-InSb 
have  always  shown  minima  near  the  resonance  points. 

In  each  of  figures  2-4,    we  show  curves  for  pure  optical  phonon 
scattering  (F=   0),    labelled  (a),    which  show  the  magnetophonon 

oscillations  most  clearly.      To  be  noted  are  the  facts  that  the  resonance 

4 
extrema  are  maxima,    as  in  the  Gurevich-Firsov  theory,      but  that  the 

maxima  lie  somewhat  to  the  low- field  side  of  the  resonance  points, 

defined  by    uu     =   uu    /N,    N  =   1,    2,    .    .      .      The  shift  is  greatest  at  low 
c  o 

temperatures  and  low  N- values.      The  pure  polar  optical  curve  at  y    =14 
also  shows  the  near- logarithmic  behavior  at  large  B,    discussed  in  the 
preceding  section. 

The  principal  result  of  adding  other  scattering  mechanisms  to 
the  polar  optical  scattering  is  to  diminish  and  finally  obliterate  the 
magnetophonon  oscillations,    as   seen  in  the  figures.      In  all  cases  we 
have  examined,    the  oscillations  have  vanished  before     T    has  reached 
unity.      One  also  sees  from  the  figures  that  there  is  no  significant 
"phase  shift,  "   due  to  these  other  mechanisms,    moving  the  minima 
toward  the  resonance  points,    which  possibility  was  considered  in  the 
introduction. 

In  figures  5  and  6,    we  exhibit  the  results  of  some  hot-electron 
calculations,    obtained  by  solving  simultaneously  the  carrier  momentum 
and  energy  steady-state  equations.     A  "  constant- current"  (actually, 
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constant  drift  velocity)  magnetoresistance    curve  is   shown  in  each 
figure.      The  ohmic  magnetoresistance  is  also  shown  for  comparison. 
The  dimensionless  drift  velocity   v   in  the  figures  is  defined  in  eq  (34). 

With  the  n-InSb  parameters  given  earlier,    figure  5   corresponds  to  a 

5 
temperature  of  20  K,    and   v=    0.  015  corresponds  to  v  =0.  11   xlO   m/sec, 

d 

or  E  =   0.  56  V/cm  at  B  =  0.      In  figure  6,    T      =  77  K,    and   v=  0.  32  corres- 

.Li 

5 
ponds  to  v  =2.  34  x  1  0   m/sec,    or  E  =   42  V/cm  at  B  =  0. 

The  most  noticeable  features  of  figures  5  and  6  are  the  sharp 

2,  3 
cusps  precisely  at  the  resonance  points.  These  are  energy-transfer 

effects   superimposed  on  the  background  of  momentum  scattering  which 

comprises  the  ohmic  curves.      Nonparabolicity  of  the  conduction  band 

and  Landau  level  broadening  would  tend  to  round  off  the  cusps.      Figure 

5,    showing  resonance  maxima,    is  interesting  in  that  it  shows  how  the 

magnetophonon  structure,    absent  in  the  ohmic  case,    can  be  brought  out 

by  application  of  nonohmic  electric  field  or  current.      This  effect  has 

38,  39  40 

been  exploited  by  Stradling  and  Wood,  and  by  Aksel'rod  et.    al. 

In  figure  5  the  scattering  is  dominated  by  ionized  impurities,    the 

-3/2 
resistivity  of  which  varies  essentially  at  T  at  both  B  =0  and  the 

quantum  limit.      Since  the  carrier  temperature  decreases  to  the  lattice 

2,  3  . 

temperature  at  the  resonance  points,  the  magnetoresistance  in- 

creases.     Figure  6  shows  resonance  cusps  which  are  minima.      At  this 
temperature,    the  optical  phonon  scattering  dominates  both  the  acoustic 
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phonon  and  impurity  scattering.      The  minima  are  the  result  of  the 
polar  optical  phonon  scattering  itself. 

Low  temperature,    hot- electron  longitudinal  magnetoresistance 

38  41  42 

in  GaAs,         CdTe,         and  InAs        do  show  maxima,    in  agreement  with 

39,  40 
figure  5.      However,    InSb  is  an  exception.  The  reason  for  this 

43 
is  not  clear,    although  Yamada  and  Kurosawa       argued  that  a  deviation 

from  the  Maxwellian  distribution  function  could  occur  under  hot- 
electron  conditions  and  could  account  for  the  minima.      Other,    physical, 

39 
arguments  can  also  be  advanced      ,    to  explain,    in  part,    the  minima. 

Further  investigation  is  in  order  on  this  point. 

In  summary,    the  displaced-Maxwellian  distribution  function 

together  with  the  simple  model  of  the  band  structure  cannot  explain  the 

observed  longitudinal  magnetoresistance  resonance  minima  which 

frequently  oceur  under  ohmic  conditions.      This  is  in  spite  of  the  fact 

that  the  B  =   0  and  quantum  limit  mobilities  largely  agree,    except  for 

numerical  factors,    whether  computed  from  the  displaced-Maxwellian, 

or  the  Boltzmann  equation  without  carrier- carrier  scattering.      It 

does  not  seem  likely  that  the  simple  band  structure  assumed  (parabolic 

and  isotropic)  is  the  principal  reason  for  this  failure,    since  apart 

from  InSb,    the  semiconductors   studied  do  not  deviate  greatly  from 

the  simple  model  (subsidiary  valleys  are  surely  not  important  at 

ohmic  conditions),    and  yet  exhibit  near-resonance  minima,    or  in 
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some  cases  an  ambiguous  situation.     More  likely,    the  lack  of  agree- 
ment lies  with  the  displaced-Maxwellian  itself.      Several  factors  point 

4 
to  this.     First  of  all,    the  original  work  of  Gurevich  and  Firsov, 

based  on  an  approximate  but  direct  solution  of  the  Boltzmann  equation, 

implies  this.      The  materials  studied  experimentally  have  usually  been 

15  -3 

fairly  pure  (  -      10        cm        impurity  concentration),    whereas  the  dis- 


placed-Maxwellian implies  a  carrier  concentration  usually  greater 

15  -3 

than  about  10        cm        (this  figure  is  quite  uncertain;  it  is  a  function 

of  the  strengths  of  all  the  carrier-phonon  and  carrier-impurity  inter- 
actions,   as  well  as  of  temperature,    and  electric  and  magnetic  field 
strengths      ).      Further,    the  work  of  Yamada  and  Kurosawa        is  a 
further  indication  of  the  inadequacy  of  the  displaced-Maxwellian  for 
the  usual  experimental  situations.      Nonohmic  magnetophonon  results 

of  the  displaced-Maxwellian,    however,    can  be  qualitatively  satisfactory, 

3 
as  noted  above  and  also  previously. 
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